Instanton size distribution: repulsion or the infrared fixed point? 
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We discuss available information about the instanton size 
distribution d(p), which comes from lattice simulations and 
the interacting instanton liquid model. Not only they are 
remarkably consistent, but one can also reproduce d(p) with 
an alternative idea, based on the infrared fixed point. We 
show, that lattice non-perturbative beta-function (for bare 
charge renormalization) also suggests such a hint. We also 
discuss whether it is possible to reconcile the instanton physics 
with reasonable large-number-of-colors limit, concluding that 
our fits for SU(2) and SU(3) theories are actually very close 
to the critical value 8tv 2 /q 2 = N c * 5.4 when it is the case. 



Instantons were discovered 20 years ago pL but only 
now their significant role in strong interaction physics 
is being widely recognized. A large set (~ 40) of 
point-to-point correlation functions was calculated in the 
simplest "instanton liquid" model M, both numerically 
H and (with certain approximations) analytically [Q. 
The results are truly spectacular: not only masses and 
coupling constants of lowest hadronic states (including 
7r, a, p, A\, N, A and others) are reproduced without ad- 
ditional parameters, but the whole correlation functions 
happen to be in agreement with those extracted from 
phenomenology || and from lattice simulations || . This 
implies that instanton-induced 't Hooft interaction ]Tl| ] 
between light quarks is dominating the interquark forces. 
Recently, glueballs were added to the list p0| , with the 
conclusion that their specific features in the instanton 
model (such as small mass and especially small size of 
the scalar glueball) can quantitatively explain available 
lattice data. 

This agreement is not accidental, as can be seen from 
the fact that the "instanton liquid" itself was "distilled" 
from lattice configurations by the "cooling" method 
|RP]j§)- In spite of the (nearly complete) loss of confine- 
ment and of the perturbative gluons, this operation does 
not affect much these correlators (and hadrons) Q . Fur- 
thermore, the main parameters (the average instantons 
size p and instanton spacing R) were found to be in- 
side 10% the same as those suggested in 1982 Q], namely 
p ~ 1/3/to, R m Ifm. 

The topic of this work is the shape of the instanton 
(plus antiinstanton)size distribution d(p), related with a 
number of important issues, both practical and theoret- 
ical ones. With all those advances at the phenomeno- 
logical front, we are still lacking answers to many major 



questions. One of them, to be discussed in this letter, 
is: Why are large — size instantons absent in the QCD 
vacuum? Alternative explanations are: (i) they are sup- 
pressed by repulsive interaction between instantons; (ii) 
the higher-order effects lead to charge renormalization 
so that their actions are large; (iii) confinement effects 
screen their gluoelectric fields. In this Letter we will go 
through this list, and will show that the first two are still 
strongly competitive. 

The first idea is the basis of the 'interacting instanton 
liquid model' (IILM), formulated as a particular statis- 
tical model amenable for numerical simulations Jl4] , p^| , 
its partition function (for N/2 instantons and N/2 anti- 
instantons) is generally given by 



Z = 
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where fli, denote the orientation, position and the size 
of pseudoparticle /, do(p) here corresponds to non — 
interacting instantons, with the gluonic interactions in- 
cluded separately. (The quark-related ones are not dis- 
cussed in this letter.) 

The simplest "hard core" model was introduced in [p_ 2| - 
Then it was shown that for the simplest "sum ansatz" |13| ] 
such a repulsion actually exists, even slightly stronger 
than it is needed. Some defects of this ansatz was cured 



in 1 14 1, where the improved trial function known as "ra- 
tio ansatz" was proposed: this reproduces phenomeno- 
logical parameters of the "instanton liquid" well. How- 
ever, numerical solution |]l6| of the "streamline" equation 
p7[ have shown that the instanton-antiinstanton valley 
leads continuously to zero fields. Thus, it seems that 
the original hopes to stabilize the ensemble at the purely 
classical level |l| are not fulfilled. Presumably quan- 
tum effects (especially subtraction of perturbative contri- 
butions, relevant for close instanton- anti-instanton pairs 
with a strongly attractive interaction) will generate the 
effective repulsion (see also ||l8| ). Meanwhile, in IILM- 
based recent studies jlO],[L9| the streamline interaction is 
supplemented by a repulsive core, with the radius fitted 
to the value of the gluon and quark condensates. 

In Fig.l (open points) we show the resulting instanton 
size distribution in IILM, for the pure gauge SU(2) (a) 
and SU(3) (b) theories. In the former case they are com- 
pared to first available lattice data ]8| (closed points). 
They are very similar, except at small sizes (instantons 
with the radius p ~ 1 in lattice units "fall through the lat- 
tice" during "cooling" ) . At large sizes both lattice data 
and "interacting liquid" give virtually identical results. 
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FIG. 1. The instanton size distribution in d=4 SU(2) (a) 
and SU(3) (b) gauge theories. The open points correspond to 
interacting instanton liquid", and the closed ones are lattice 
results By Michael and Spenser, lQ 4 ,A/g 2 = 2.4 for squares 
and 24 4 , 4/g 2 = 2.5 for dots (the dotted and dashed lines just 
guide the eye). Solid curves correspond to the parametriza- 
tion discussed in the text. Units are in "femtometers" , de- 
fined for lattice data by the scalar glueball mass defined as 
m + = 1.7"GeV"; and by 1/Apv for the instanton model. 

It is well known that (provided the action is large 
enough S e ff » 1) one can use the semiclassical the- 
ory and derive the following generic expression 



d(p) 



2d 

0(^+1; 



Sl N ™M(p)exp[-S(j>)] 



(2) 



where 2 accounts for instantons and anti-instantons, d is 
the space-time dimension, and N zm is the number of zero 
modes (4iV c for Yang-Mills fields). The value of the con- 
stant (which depends on renormalization method) was 
determined for the Yang-Mills instantons in the classical 
work by 't Hooft [ pJ| . In order to make our discussion 
below simpler, we use C/ = 0.466/[(7V c - l)!(iV c - 2)!], 
absorbing the N c dependent factor into a new Lambda 



Multiple arguments in favor of existence of the infrared 
fixed point (or "freezing" of the coupling constant) have 
been many times made before, based on a variety of 
phenomenological models (one may find a list of refer- 
ences e.g. in a recent paper pQ], which came to the 
same conclusion via analysis of higher-loop corrections 
to <r(e + e~ —> hadrons) ). If so, at large p the action is 
p-independent, and the following simple prediction fol- 
lows: d(rho) ~ p~ d+1 , where d is just the space-time 
dimension. Are available data in agreement with it? 

In terms of statistical accuracy and the widest range 
studied, the best lattice measurements are those for d=2 
0(3) sigma model |2l|] . For large sizes the result is d(p) ~ 
P 3 @j m perfect agreement with the idea of "frozen" 
S(p). 

Furthermore, for the d=4 SU(2) gauge theory (see 
Fig. 1(a)) we have fitted the action using the following 
simple parametrization, a standard two-loop expression 
for the charge (where b = ^-N c ,bi = ^-N^) 



8tt 2 

Tip) 



= bgL + b\logL 



with a "regularized" log 



L = -log[{\y + C] 
P pA 



(3) 



(4) 



It has two new parameters, C and p, describing where and 
how rapidly the "freezing" occurs. Solid line in Fig. la 
shows such a fit, with A inst =.66 /to -1 , p=3.5, C=4.8. 
For the SU(3) gauge theory the data for d(p) itself are still 
missing, but the total instanton density and average size 
were measured in 0: J dpd(p) « 1.3/m -4 ,/? = 0.35/to. 
One may fix the parameters to reproduce those two inte- 
grals: the corresponding solid curve is shown in Fig. 1(b). 
(Parameters in this case are A,- ns j = .70/m , p=3.5, 
C=5.0. ). Although it is not a fit to open points (the 
interacting liquid) and is based on different physics, the 
curve agrees with the points perfectly. Now we plot both 
fits it in terms of the action (Fig. 2a) one can see that 
it indeed means rather rapid "freezing" of the coupling 
constants. We do not know why such freezing may occur: 
but it is remarkable that it should happen where the ac- 
tion itself is large, S ~ 10 — 20 (so that one should not 
question the semiclassical theory). 
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FIG. 2. The upper part shows the instanton action S(p) 
versus p (in "fm") according to the parametrization used in 
Fig.l. The lower part compare the fitted values of the "fixed 
point" actions S* (divided by the number of colors N c = 2, 3) 
with the critical point in the large N c limit (star). 

Let us now comment on the alternative (iii) mentioned 
above, namely that the cutoff of the large size instantons 
is related with confinement. Two arguments against its 
dominance can be given. Although no specific formulae 
are known (they should depend on the unknown confine- 
ment mechanism), the confinement- related correction to 
the action should have the form 5S con f(Kp 2 ) where K is 
the string tension. Most probably it is a regular function, 
expandable in powers of its argument: this leads to an 
exponential cutoff at large p rather than power-like one, 
as suggested by Fig.l. The second: studies of the in- 
stanton sizes at non-zero temperatures p2| have recently 
been made. Near the deconfinement transition (where 
K(T) — > 0) the average instanton size does not grow, 
while one should expect it to happen (as 1/K 1 / 2 ) if con- 
finement is involved. (One may investigate this question 
further by a number of methods: in particular, by check- 
ing whether the so called "abelian projected monopole 
loops" |2(j] are or are not correlated with instantons in 
the vacuum ensemble.) 

Now, is the existence of the infrared fixed point consis- 
tent with non-perturbative beta functions used in the lat- 
tice studies? Recall that the bare coupling g(a) is fixed at 
the input of simulations, then one calculates observables 



(hadronic masses, etc) and fixes the physical magnitude 
of the lattice scale a. Reversing the function, one gets 
the renormalized bare charge g(a), subject for (i) univer- 
sality test (its independence on the particular observable 
used) and (ii) comparison with the expected perturba- 
tive behaviour ("asymptotic scaling"). The non-trivial 
fact found (see e.g. (2^j2^|) is that (i) extends beyond 
(ii). 

A sample of lattice data for the SU(3) lattice gauge 
theory, without quarks p4| ] (open points) and with two 
massless quark flavors (closed points) j25j shown in Fig. 3. 
Those are usually presented in form of the derivative of 
g(a), the famous Gell-Mann-Low beta function. In order 
to show various theories, one can normalize the derivative 
to its asymptotic (g — > 0) perturbative value: thus we 
plot the following ratio 



M9) = 



rf(l/g 2 ) rf(l/g 2 ) 
dlog(a) l/[ dlog{a) 1 



(5) 



which tends to 1 the right hand side of Fig. 3. As one 
penetrates into the non-perturbative region the data dis- 
play a deep drop of Rp at 6/g w 6, with subsequent 
turn upward. Both are very rapid: the latter turn is 
known in literature as a transition to a "strong coupling" 
regime. This drop may be an indication for the infrared 
fixed point, while the upturn is definitely a lattice arti- 
fact. (Rapid turnovers are common when a renormaliza- 
tion trajectory is going toward the fixed point, and then 
misses it closely.) We conjecture that for lattice actions 
which are better than Wilson one, the corresponding beta 
function would drop further. 

However, even with standard Wilson actions one may 
investigate other definitions of the renormalized charge, 
namely those related with the renormalized background 
fields. This method is quite popular in the perturbative 
context p^] but (to our knowledge) has not yet been used 
on the lattice. It is clear which field is the best to try. 
First, to avoid complications with external current, one 
should better take a "self-supporting" classical field, such 
as D^GL, = 0. Second, topology adds additional stabil- 
ity: thus instanton is the most natural classical back- 
ground. Third, its renormalized action depend on one 
parameter p, S e ff — 87r 2 / ' g 1 \p) , and thus this expres- 
sion can be used to study the non-perturbative charge 
renormalization. 
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FIG. 3. The non-pertubative beta function for the SU(3) 
lattice theory with Wilson action. Open points are taken 
from R.Gupta review (which is an analysis of the data set 
from Lapage and Mackenzie), the solid points are from Blum 
et al. The solid line is a fit discussed in the text. 

One can put an instanton on the lattice and then "heat 
it up" , performing standard updates: such studies have 
been made by DiGiacomo et al [p3[ , establishing renor- 
malization of the topological susceptibility. However, for 
the proposed goal is not enough to keep a topological 
charge (as done in ref. one has to preserve the 

chosen value of p. This can be achieved in two ways: (i) 
while updating a link, one may keep the quantum field a M 
orthogonal to the dilatational zero mode SA^/Sp p9| ; or 
(ii) one may use a modified lattice action containing the 
two-plaquette operators with parameters tuned to make 
any given size the classical minimum of the lattice action. 
Either way, the main problem is to get high statistics 
measurements of the effective action, after subtraction of 
the usual "average plaquette" is made. 

The next topic addressed in this Letter is the fate of 
instantons in the large N c limit. It was first addressed by 
Witten J27J, who has formulated the following dilemma: 
either (i) instantons are not dynamically relevant, while 
all observables have regular 1 /N c expansion as suggested 
by perturbative expansion, or (ii) instantons play some 
role in the real world, but are exponentially suppressed at 
large N c . Back in 1979, Witten has argued in favor of (i), 
using analogies with some d=2 models, but today there 
is no doubt about significance of the instanton-induced 
effects at N c = 3. However, it is hard to accept (ii) also, 
because then all the 1/N C development is undermined. 

This dilemma still allows for one loophole Q , an alter- 
native (iii), which allows to reconcile instanton physics 
with smooth large N c limit. Let us look at its conse- 
quences and compare them with the empirical trends 
found for N c — 2, 3. The semiclassical formula (2) in 



the large N c limit contains factorial terms in the de- 
nominator. The action should grow linearly with N c : 
S(p) = N c s(p) and therefore these factorials are exactly 
cancelled by the gauge zero modes. The next level are the 
exponential terms, which with our normalization look as 
follows 

d(p) ~ exp{N c [2 - s(p) + 2logs(p)}} (6) 

The limit clearly depends on the sign of the bracket. If 
there exist the fixed point s(p) — > s* , as we advocated 
above, it is important whether it is larger or smaller than 
the critical value s* = 5.4, the root of the bracket. 

Let us compare the values of the fixed actions for 
N c = 2,3 obtained from the fits above to the critical 
value (shown by the star) in Fig. 2(b). The first obser- 
vation is that for two fitted cases the limiting action S* 
indeed grow linearly with N c . Second, the constant is 
somewhat larger than the critical value, although this 
is so to say inside the error bars (we have assign those 
to the points, which relate to the uncertainty due to in- 
stanton interactions) . In general, it seems quite plausible 
that at N c — > oo instantons will take care of themselves 
and their density and role will remain roughly constant. 

But should we expect the limit N c — > oo to be smooth? 
The answer is negative (l|: at larger N c an "instanton 
liquid" should become solid, because a growing action 
(which depends on relative positions) is analogous to a 
decreasing temperature. Thus, one should expect that 
both color and translational symmetries would be spon- 
taneously broken. In |l5[ ] such phase transition was in 
fact found for N c = 3, but at non — physical instanton 
density (about 60 times the physical one). For larger N c 
the work is in progress. Lattice studies of theories with 
N c > 3 would certainly be of great value. 

Finally, let us mention the practical aspect of the stud- 
ies of d(p): its measurements at small p < 0.2 fm is 
potentially the source of by far the most accurate mea- 
surements of Aqcd- As it is well known, here d(p) ~ 
^QCDP^ b °~ d ~^ : tnus even relatively poor accuracy (say 
50%, compare to dots at Fig. la) leads to Aqcd with ac- 
curacy of about 5%, much better than measured today. 

In summary, we have argued that existing data for the 
instanton size distribution can be equally well be de- 
scribed both by the interacting instanton liquid model 
(IILM) and by the existence of the infrared fixed point. 
If the second is the case, it seems plausible that this fixed 
point is not far from the critical point, separating sce- 
narios in which instantons are suppressed or grow in the 
N c — > oo limit. Better measurements of d{p) on the lat- 
tice (especially at N c = 3 and larger) are badly needed. 
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